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Abstract 

Many-body fermion systems are important in many branches of physics, including condensed 
matter, nuclear, and now cold atom physics. In many cases, the interactions between fermions 
can be approximated by a contact interaction. A recent theoretical advance in the study of these 
systems is the derivation of a number of exact universal relations that are predicted to be valid for 



all interaction strengths, temperatures, and spin compositions [iHSl]. These equations, referred to 
as the Tan relations, relate a microscopic quantity, namely, the amplitude of the high-momentum 
tail of the fermion momentum distribution, to the thermodynamics of the many-body system. In 
this work, we provide experimental verification of the Tan relations in a strongly interacting gas of 
fermionic atoms. Specifically, we measure the fermion momentum distribution using two different 
techniques, as well as the rf excitation spectrum and determine the effect of interactions on these 
microscopic probes. We then measure the potential energy and release energy of the trapped gas 
and test the predicted universal relations. 
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In 2008, Shina Tan derived a number of universal relations for an interacting Fermi gas 
with short-range, or contact, interactions iHSj. These relations are very powerful because 
they connect microscopic quantities, such as the momentum distribution of the fermions, to 
macroscopic quantities, such as the total energy of the system. Furthermore, the relations are 
universal in that they do not depend on the details of the interparticle potential that gives rise 
to the interaction, nor do they depend on the state of the system, which could be an exotic 
Fermi superfiuid, a normal Fermi liquid, or even a simple two-body state such as a diatomic 
molecule composed of two fermionic atoms. At the heart of the universal relations is a single 
quantity, which Tan termed the contact. The contact is defined as the amplitude of the high- 
k tail of the momentum distribution n{k), which was previously predicted to scale as 
for an interacting Fermi gas j?]. Remarkably, it can be shown that the contact encapsulates 
all of the many-body physics Therefore, it is predicted that by measuring this tail 
of the momentum distribution and then applying universal relations, one could determine 
many other properties of the system. Two recent papers report the contact for a strongly 
interacting Fermi gas^^ extracted from photoassociation measurements jsl, O] and inelastic 
Bragg spectroscopy [9] and compare the results with theoretical predictions for the BCS- 
BEC crossover. Here, we present a series of measurements that not only measure the contact 
in the BCS-BEC crossover with two different techniques, but moreover test the Tan relations 
experimentally by comparing measurements of both microscopic and macroscopic quantities 
in the same system. Our results directly verify the universal relations by exploiting the fact 
that while the value of the contact depends on the many-body state and on parameters such 
as temperature, number density, and interaction strength, the universal relations do not. 

Our measurements are done in an ultra cold gas of fermionic atoms confined in a 
harmonic trapping potential. We cool the gas to quantum degeneracy in a far-detuned 
optical dipole trap as described in previous work 10|. The trap is axially symmetric and 
parameterized by a radial trap frequency, which varies for these data from Ur = 2tt ■ 230 
to 27r ■ 260 Hz and an axial trap frequency, which varies from = 27r ■ 17 to 27r • 21 Hz. 
We obtain a 50/50 mixture of atoms in two spin states, namely the \ f,mf) = |9/2, —9/2) 
and 1 9/2, —7/2) states, where / is the total atomic spin and is the projection along the 
magnetic- field axis. Our final stage of evaporation occurs at a magnetic field of 203.5 G, 
where the s-wave scattering length, a, that characterizes the interactions between atoms in 
the |9/2, —9/2) and |9/2, —7/2) states is approximately 800 Oq, where Qq is the Bohr radius. 
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At the end of the evaporation, we have 10^ atoms per spin state at a normahzed temperature 
^ = 0.11 ± 0.02 where the Fermi temperature corresponds to the Fermi Energy, Ep = 
kbTp = hu}{6NY^^ where N is the total atom number in one spin state and to = {lo'^uJzY^^ and 
kb is the the Boltzmann constant. After the evaporation we increase the interaction strength 
adiabatically with a slow magnetic-field ramp to a Fano-Feshbach scattering resonance. 

The momentum distribution of the fermions, n{k), is predicted to scale as 1/k^ at high 
k, with the contact being the coefficient of this high momentum tail. Following Tan [2^, we 
define the integrated contact per particle for the trapped gas, which we will refer to simply 
as the contact, using 

C = lim k^n{k). (1) 

fc— >oo 

Here, k is the wave number in units of the Fermi wave number, kp = , and n{k) 
for a 50/50 spin mixture is normalized such that /q°° -^^d^k = 0.5. Note that the contact 
is expected to be the same for both spin states in the interacting Fermi gas, even in the 
case of an imbalanced spin mixture. Theoretically, the contact is defined in the limit of 
1 << k « l/(fci?ro), where tq is the range of the interatomic potential. Using a typical 
value of kp = 2^%^ for our trapped gas and the van der Waals length, tq = 60ao, we 
find l/ikpro) = 37. 

We directly measure n{k) using ballistic expansion of the trapped gas, where we turn off 
the interactions for the expansion. We accomplish this by rapidly sweeping the magnetic 
field to 209.2 G where a vanishes, and then immediately turning off the external trapping 
potential HI]. We let the gas expand for 6 ms before taking an absorption image of the 
cloud. The probe light for the imaging propagates along the axial direction of the trap and 
thus we measure the radial momentum distribution. Assuming the momentum distribution 
is spherically symmetric, we obtain the full momentum distribution with an inverse Abel 
transform. 

Fig. [iK shows an example n{k) for a strongly interacting gas measured with this technique. 
For this data, the dimensionless interactions strength (kpa)'^, is —0.08 ±0.10. Empirically, 
we find that the measured n{k) exhibits a, 1/k^ tail and we extract the contact C from 
the average value of k'^n{k) for k > kc where we use kc = 1.85 for {kpo}'^ > —0.5, and 
kc = 1.55 for {kpa)^^ < —0.5. One issue is whether or not the interactions are switched 
off sufficiently quickly to accurately measure the high-A; part of the n{k). The data in 
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FIG. 1: Extracting the contact from the momentum distribution and rf Uneshape. (a) 

Measured momentum distribution for a Fermi gas at = —0.08 it 0.10. Here, the wave number 
k is given in units of kp, and we plot the normalized n{k) multiplied by k^. The dashed line 
corresponds to C = 2.15 ± 0.2, which is obtained by averaging k'^n{k) for k > 1.85. (Inset) The 
measured value for C depends on the rate of the magnetic-field ramp that turns off the interactions 
before time-of-flight expansion, (b) RF lineshape measured for a Fermi gas at = —0.03 it 0.10. 
Here, is the rf detuning from the single-particle Zeeman resonance, given in units of Ep /h. We 
plot the normalized rf lineshape multiplied by 2^/^7r^z^^/^, which is predicted to asymptote to C 
for large i'. Here, we obtain C = 2.13 it 0.2 (dashed line) from an average of the data for > 5. 

Fig. [1^ were taken using a magnetic-field sweep rate of approximately 1.4^ to turn off 
the interactions for the expansion. In the inset to Fig. [1^, we show the dependence of the 
measured C on the magnetic-field sweep rate. Using an empirical exponential fit (line in Fig. 
[It inset), we estimate that at our typical sweep rate of approximately 1.4^ our measured 
C is systematically low by about 10%. We have therefore scaled the contact measured with 
this method by 1.1. 

The contact is also manifest in rf spectroscopy, where one applies a pulsed rf field and 
counts the number of atoms that are transferred from one of the two original spin states 
into a third, previously unoccupied, spin state We transfer atoms from the |9/2, —7/2) 
state to the |9/2,— 5/2) state. It is predicted that the number of atoms transferred as a 
function of the rf frequency, scales as z/~'^/^ for large u, and that the amplitude of this 



high frequency tail is 
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4l5l|. Here, = is the single-particle spin-fiip resonance. 



and V is given in units of Ep / h. This prediction requires that atoms transferred to the third 
spin-state have only weak interactions with the other atoms, so that "final state effects" are 
negligible 15|-|22|. as is the case for '^^K atoms. In Fig. [T]d, we plot a measured rf spectrum 
multiplied by 2^/^7r^z/'^/^. The rf spectrum, r(z/), is normalized so that the integral over the 
rf lineshape equals 0.5. Empirically, we observe the predicted l/v^/"^ behavior for v > uc- To 
obtain the contact we average 2^/^7r^z/^/^r(z/) for u > uc where z/c = 5 for [kpo)^^ > —0.5, 
and i^c = 3 for {kpa)^^ < —0.5. 

The connection between the tail of the rf spectrum and the high-fc tail of the momentum 
distribution can be seen in the Fermi spectral function, which can be probed using photoe- 
mission spectroscopy for ultra cold atoms |p^]. Recent photoemission spectroscopy results 



on a strongly interacting Fermi gas 



23| revealed a weak, negatively dispersing feature at 



high k that persists to temperatures well above Tp. This feature was attributed to the effect 
of interactions, or the contact, consistent with a recent prediction that the l/k"^ tail in n{k) 
should correspond to a high- A; part of the spectral function that disperses as —k'^ [24]. Atom 
photoemission spectroscopy, which is based upon momentum-resolved rf spectroscopy, also 
provides a method for measuring n{k). By integrating over the energy axis, or equivalently, 
summing data taken for different rf frequencies, we obtain n{k). This alternative method 
for measuring n{k) yields results similar to the ballistic expansion technique, but avoids the 
issue of magnetic-field ramp rates. 

In Fig. [2] we show the measured contact for different values of the dimensionless inter- 
action strength, 1/kpa. Here, the contact is extracted using the three different techniques 
described above to probe two distinct microscopic quantities, namely the momentum dis- 
tribution and the rf lineshape. We find that the amplitude of the l//c^ tail of n{k) and the 
coefficient of the tail of the rf spectra yield consistent values for C. The solid line is 

a prediction for the contact that was reported in Fig. 1 of Ref . \^ . This prediction consists 
of the BCS limit, interpolation of Monte Carlo data near unitarity, and the EEC limit for a 
trapped gas at zero temperature and uses a local density approximation. 

Remarkably, the Tan relations predict that the contact, as revealed in probes of the 
microscopic behavior of the gas, is directly connected to the thermodynamics of the gas. To 
test the Tan relations, we measure the potential energy and release energy of the cloud. The 
total energy of the trapped gas divided by the number of particles, E, is the sum of three 
contributions, the kinetic energy T, the external potential energy V, and the interaction 
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FIG. 2: The contact. We measure the contact, C, as a function of (kpa)'^ using three different 
methods. Filled circles correspond to direct measurements of the fermion momentum distribution 
n{k) using a fast magnetic-field sweep to project the many-body state onto a non-interacting state. 
The momentum distribution can then be measured in ballistic expansion. Open circles correspond 
to n{k) obtained using atom photoemission spectroscopy measurements. Stars correspond to the 
contact obtained from rf spectroscopy. The values obtained with these different methods show 
good agreement. The contact is nearly zero for a weakly interacting Fermi gas with attractive 
interactions (left hand side of plot) and then increases as the interaction strength increases to the 
unitarity regime where (kpa)^^ = 0. The line is a theory curve obtained from Ref. 



energy /. Because the momentum distribution has a tail, the kinetic energy obtained 
by integrating over the momentum distribution diverges. However, the release energy T + 
/, which can be measured in the usual time-of-flight expansion without turning off the 
interactions, is not divergent. 

We measure V by imaging the spatial distribution of the atom cloud, similar to measure- 



ments that were performed in Ref. 



251]. We allow the cloud to expand for 1.6 ms to lower 



the optical density and then image along one of the radial directions in order to see the 
density distribution in the axial direction. Because the expansion time is short compared 
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to the axial trap period (40 times shorter) the density distribution in the axial direction 
reflects the in-trap axial density distribution. The potential energy per particle, in units of 
Ep, is then V = jj^^mco'^{z'^), where (z^) is the mean squared width of the cloud in the 
axial direction and we have assumed that the total potential energy is distributed equally 
over the three axes. 

To measure the release energy T + / we turn off the trap suddenly and let the cloud 
expand for t = 16 ms (with interactions) before imaging along one of the radial directions; 



this is similar to measurements reported in Ref. [26|]. The total release energy is the 
sum of the release energy in the two radial directions and the release energy in the axial 
direction. For the radial direction, the release energy per particle, in units of Ep, is simply 
Tr + Ir = ]v^l'^'^ where t is the expansion time and (y^) is the mean squared width of the 
expanded cloud in the radial direction. For the axial direction, the expansion is slower and 
the expanded cloud may not be much larger than the in-trap density distribution. This is 
especially true near the Feshbach resonance where the cloud expands hydro dynamically 
Accounting for this, the axial release energy is + = -jy^^m ^^ ^2 ^° ; where Zq is mean 
squared axial width of the in-trap density distribution. We extract the mean squared cloud 
widths from surface fits to the images, where we fit to a finite temperature Fermi Dirac 
distribution while minimizing the difference in energy between the raw data and the fit. 
Rather than being theoretically motivated, we simply find empirically that this functional 
form fits well to our images. To eliminate systematic error due to uncertainty in the trap 
frequencies and imaging magnification, we measure the release energy and potential energy 
of a very weakly interacting Fermi gas at ^ = 0.11, where T + I and V for an ideal Fermi 
gas is OAOEp. We then use the ratio of OAOEp to our measured values as a multiplicative 
correction factor that we apply to all of the data. This correction is within 5% of unity. 
For the point with > we take into account the binding energy of the molecules in our 
tabulation of the release energy T + / by adding —l/lkpa)"^. We show our data for the V 
and T + I versus {kpa)~^ in the inset of Fig. [3l 

We can now test the predicted universal relations connecting the l/fc^ tail of the momen- 
tum distribution with the thermodynamics of the trapped Fermi gas. We first consider the 
adiabatic sweep theorem [l|, 

which relates the contact C to the change in the total energy of the system when the 
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FIG. 3: Testing the adiabatic sweep theorem. (Inset) The measured potential energy, V, and 
release energy, T + I, per particle in units of Ep is shown as a function of l/kpa. (Main) Taking 
a discrete derivative of the data shown in the inset, we find that 27r ^^_-^^^^^^^ (•) agrees well with 
the contact C measured from the high-Zc tail momentum distribution (o). 

interaction strength is changed adiabatically. To obtain the energy per particle, E, we sum 
the values for T + / and V shown in the inset of Fig. [31 To test the adiabatic sweep theorem, 
we find the derivative, , simply by calculating the slope for pairs of neighboring 

points in the inset to Fig. |3l In the main part of Fig. HI we compare this point-by-point 
derivative, multiplied by 27r, to C obtained from the average values of the data shown in 
Fig. |2](o). Comparing these measurements of the left and right sides of Eqn. [21 we find good 
agreement and thus verify the adiabatic sweep theorem for our strongly interacting Fermi 
gas. 

A second universal relation that we can directly test is the generalized virial theorem (2|, 

C 



E-2V = T + I-V 



(3) 



which relates the difference between the release energy and the potential energy to the 
contact. Eqn. [3] is predicted to be valid for all values of the interaction strength {kpa)'^- 
This generalized virial theorem reduces to E — IV = for the ideal gas, where / = 0, as well 
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FIG. 4: Testing the generalized virial theorem. The difference between the measured release 
energy and potential energy per particle T+I — V is shown as open circles. This corresponds to the 
left hand side of Eq. [3j Filled circles show the right hand side of Eq. [3] obtained from the average 
values of the contact shown in Fig. [21 The two quantities are equal to within the measurement 
uncertainty, which is on order of O.OlE'i?. 

as for the unitarity gas, where {kpo}^^ = 0. This result for the unitarity gas was previously 



verified in Ref. 
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|. Here, we test Eqn. |3]for a range of interaction strengths. In Fig. IHwe 
plot the measured difference T + I — V versus {kpa)~^ along with ^J^^^ , where we use our 
direct measurements of C. We find that these independent measurements of the left and 
the right sides of Eqn. |3] agree to within our error, which is roughly 1% of the Fermi energy. 
It is interesting to note that the measured energy difference T + / — is small (in units of 
Ep), so that even a Fermi gas with a strongly attractive contact interaction nearly obeys 
the non-interacting virial equation. 

In conclusion, we have measured the integrated contact for a strongly interacting Fermi 
gas and demonstrated the connection between the tail of the momentum distribution 
and the high frequency tail of rf spectra. Combining a measurement of C vs (kpa)'^ with 
measurements of the potential energy and the release energy of the trapped gas, we verify 
two universal relationships 2l, ISj , namely the adiabatic sweep theorem and the generalized 
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virial theorem. These universal relations, which have now been experimentally confirmed, 
represent a significant advance in the understanding of many-body quantum systems with 
strong short-range interactions. Furthermore, these connections between microscopic and 
macroscopic quantities could be exploited to develop novel experimental probes of the many- 
body physics of strongly interacting quantum gases. 
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